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1.  Introduction 

Rough-wall  turbulent  flows  are  more  common  in  engineering  application 
than  smooth-wall  turbulent  flows.  Modification  of  mean  flow  and  turbu- 
lence property  have  been  established  by  enormous  experiments.  However, 
details  of  mechanism  of  rough-wall  turbulence  has  been  understood  little 
because  spatial  resolution  is  limited  in  experiments.  Meanwhile,  Direct  Nu- 
merical Simulation  (DNS)  of  a rough-wall  turbulent  flows  which  is  capable 
of  giving  high  resolution  data  requires  prohibitively  heavy  computer  power 
and  accordingly,  no  other  DNS  data  are  available  than  those  published  by 
the  present  authors  (Miyake  et  al,  1999).  Our  first  DNS  considered  sand- 
grain  roughness  whose  effect  was  implemented  by  profile  drag  based  on 
Stokes  drag  and  could  successfully  reproduce  experimentally  established 
rough-wall  turbulent  flow  such  as  downward  shift  of  straight  line  of  loga- 
rithmic mean  velocity  distribution  and  vanishing  of  viscous  sublayer.  It  was 
confirmed  that  the  layer  adjacent  to  the  wall  up  to  several  tens  in  wall  unit 
where  smooth-wall  turbulence  exhibits  autonomous  property  independent 
on  the  layer  above  it,  is  taken  over  by  the  layer  having  property  of  log- 
arithmic layer,  in  rough-wall  turbulence.  While  quasi-streamwise  vortices 
play  major  role  to  generate  high  turbulent  shear  stress  in  this  near-wall 
layer  in  smooth-wall  turbulent  flow,  roughness  destroy  this  vortical  sys- 
tem and  consequently,  different  mixing  system  which  replaces  the  role  of 
quasi-streamwise  vortices  should  be  found  in  rough- wall  layer.  Present  work 
intends  to  investigate  the  turbulent  mixing  in  the  layer  close  to  the  wall  of 
rough- wall  turbulence  by  a DNS  of  more  sound  numerical  conditions,  i.e., 
without  using  any  model  for  roughness  element. 
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Table  1.  Computational  conditions 


Streamwise 

Wall  normal 

Spanwise 

Length 

2t xH 

2 H 

7 tH 

(wall  unit) 

(942) 

(300) 

(471) 

Grid  number 

256 

257 

128 

Reynolds  number 

Rer 

= Hut/v  = 150 

2.  Computational  procedure 

2.1.  CHANNELS 

Since  implement  of  a roughness  element  in  the  calculation  requires  enor- 
mously high  computational  power,  only  one  wall  is  assumed  rough  and  the 
other,  smooth.  The  configuration  of  the  channel  is  illustrated  in  Fig.l.  A 
Cartesian  coordinate  system  is  employed,  in  which  x,  y,  z are  streamwise, 
wall-normal,  and  spanwise  directions,  respectively,  and  the  velocity  compo- 
nents in  the  respective  directions  are  denoted  by  u,  v,  w. 

Computational  conditions  such  as  size  of  computational  domain,  grid 
number  and  Reynolds  number  are  given  in  Table  1.  The  global  Reynolds 
number  is  defined  as  Rer  = uTH/v  where  uT  is  global  mean  friction  velocity 
and  is  calculated  from  pressure  gradient  dp/dx  by  uT  = \(H/p)(dp/dx)]1^2, 
H is  half  channel  width  and  v and  p are  kjnematic  viscosity  and  density, 
respectively.  In  the  following,  superscript ! !)  denotes  quantities  averaged 
in  a plane  parallel  to  the  wall  ( x-z  plane),  superscript  ('),  fluctuation  from 
this  mean  and  superscript  (+),  non-dimensionalized  quantities  by  uT  which 
is  mean  friction  velocity  and  v.  The  spatial  discretization  is  a Fourier  series 
expansion  in  x and  z for  the  sake  of  periodicity  and  a Chebyshev  polynomial 
expansion  in  the  normal  direction.  Time  advancement  is  done  by  a semi- 
implicit  method  ( Crank-Nicolson  scheme  for  viscous  term  and  third  order 
Adams-Bashforth  scheme  for  the  convective  term). 

2.2.  ROUGHNESS  ELEMENT 

The  roughness  is  composed  of  periodically  arranged  two-dimensional  span- 
wise  ribs  whose  cross  section  is  a square  of  a side  hf  = 20  and  separated  in 
streamwise  direction  by  Ax+  = 118.  The  boundary  condition  on  the  solid 
wall  is  treated  rigorously  to  satisfy  non-slip  flow.  Namely,  at  each  grid  point 
xs  on  every  rib  surface,  external  force  fv  is  added  so  as  to  make  the  flow 
vanishing  there.  The  magnitude  of  fv  is  determined  iteratively  to  satisfy 
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( ):Grid  Number 


Figure  1.  Configuration  of  a 
channel  with  a rough  wall  of  2- 
dimensional  ribs. 


Figure  2.  Distribution  of  mean  veloc- 
ity u+  normalized  by  mean  friction 
velocity  uT<r  on  the  rough  wall  and 
v. 


fv(xs,t)  = av  f v(xs,t')dt'  + (3vv(xs,t)  (1) 

Jo 

where  av  and  (3V  are  negative  constants  to  prevent  the  computation  from 
instability  and  in  this  calculation,  they  are  av  = — 4x  105  and  /3V  = — 6xl02. 
These  procedures  are  proposed  by  Goldstein  et  al.  (1995). 

3.  Results 

3.1.  MEAN  FLOW  PROPERTY 

Mean  velocity  distributions  are  given  in  Fig.2.  In  order  to  make  the  straight 
line  of  logarithmic  part  parallel  to  the  one  established  for  the  ordinary 
smooth  wall,  correction  of  the  location  of  wall  surface  is  required.  In  this 
case,  upward-shifting  by  A y+  = 7.0  is  needed.  Magnitude  of  the  downward 
shift  of  the  straight  lines  from  the  case  of  the  ordinary  smooth  wall  is  within 
the  scatter  of  experimental  data  (Raupach  et  al.,  1991).  Corresponding 
shear  stress  distribution  is  given  in  Fig.3.  The  height  Hr  from  the  rough- 
wall  where  total  shear  stress  vanishes  is  defined  hereafter  as  the  equivalent 
half  width  of  the  rough- wall  channel.  The  width  Hr  is  larger  than  H because 
total  drag  on  the  rough  wall  is  larger  than  on  the  smooth  wall.  The  defect 
of  total  shear  stress  in  the  layer  close  to  the  rough  wall  found  in  Fig.3  is 
covered  by  pressure  drag  of  ribs.  It  is  noted  that  local  high  viscous  shear 
stress  fwr/ p = u(du/dy ) is  observed  in  the  layer  near  the  top  of  the  rib 
surface. 

3.2.  MODIFICATION  OF  LARGE-SCALE  FLOW 

The  modification  of  large-scale  flow  by  rough-wall  is  considered  in  this  sec- 
tion, since  it  is  substantial  for  global  mixing  . Figure  4 shows  the  ratio  of 
turbulent  shear  stress  —u'v1  to  turbulent  kinetic  energy  k,  r = —u’v1  /k. 
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Figure  3.  Distribution  of  shear 
stress  across  channel,  normalized 
by  overall  mean  friction  velocity 
uT  and  v. 


Figure  4.  Distribution  of  structure 
parameter  r = —u'v'/k. 


0 50  100  150  200 


(b)  2D  rough  wall 


Figure  5.  Iso-contour  of  spanwise  correlation  factor  of  streamwise  velocity  fluctu- 
ation. 


It  is  known  to  be  around  0.27  in  the  logarithmic  layer  and  the  larger 
the  Reynolds  number,  the  thicker  the  constant  r-layer.  Present  simulation 
shows  that  the  wall  roughness  makes  this  layer  to  reach  closer  to  the  wall, 
as  expected. 

As  mentioned  earlier,  the  profile  drag  is  dominant  on  the  rough  wall. 
However,  in  the  region  near  the  top  of  the  roughness  elements,  the  shear 
stress  must  be  sustained  mostly  by  turbulent  shear  flows.  Therefore  , strong 
turbulent  mixing  is  required  there.  While  quasi-streamwise  vortices  play 
dominant  role  for  strong  mixing  near  the  smooth  wall,  rough-wall  is  ex- 
pected to  have  different  mechanism  for  the  strong  mixing.  Figures  5 show 
the  contour  of  spanwise  correlation  of  streamwise  velocity  fluctuation.  In 
place  of  streaks  of  mean  spanwise  separation  A z+  ~ 100  which  is  a dis- 
tinguished feature  of  the  buffer  layer,  large  scale  streaks  of  much  larger 
mean  spanwise  separation  than  in  the  buffer  layer  appears.  These  large- 
scale  streaks  are  characteristic  of  wall  turbulence  of  a large  Reynolds  num- 
ber and  accordingly,  of  the  layer  close  to  the  wall  near  the  rough  surface. 
These  observation  suggests  that  the  turbulence  property  normally  found  in 
the  layer  away  from  the  wall  prevails  closer  towards  the  wall  in  a rough-wall 
case.  Namely  it  is  suggested  that  a large-scale  mixing  takes  place  actively 
over  the  roughness. 
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(a)  Original  y* 


Figure  7.  Distribution  of  separately 
Figure  6.  Iso-  surface  of  the  second  averaged  turbulent  shear  stress  due 
invariance  of  velocity  gradient  tensor  40  larSe  and  small-scale  flow. 

Q'  (Q'+  > 1.8  x 10“4). 


Fig.  8 Extracted  vortex  structure  LSE  with  the  detector  in  the  plane  at  y+  = 25. 
(a)  detector  at  near  the  rib,  (b)  position  of  detector  unconditioned.  ( Q,+  > 1.8  x 
1(T4). 

3.3.  TURBULENT  MIXING  DUE  TO  COHERENT  VORTICES 

To  investigate  the  contribution  of  coherent  vortices  of  different  scale  to  tur- 
bulent shear  stress,  a filter  is  used  to  decompose  the  flow.  The  decomposing 
filter  is  of  wavelengths  of  <5+  « 400, 6 + sa  100  in  x,  z direction,  respectively. 
Figures  6 compares  (a)  whole  scale  flow  and  (b)  large-scale  eddies  discarded 
flow,  by  the  iso-surfaces  of  the  second  invariance  of  velocity  gradient  tensor 
Q'+  = 1.8  x 10-4.  It  is  observed  that  the  coherent  flow  mainly  consists  of 
small-scale  eddies  of  above-defined  length.  Figure  7 which  shows  the  contri- 
bution of  each  component  to  turbulent  shear  stress  demonstrates  that  the 
shear  stress  is  composed  dominantly  of  the  small  scale  vortices  in  near-wall 
region  of  both  smooth-  and  rough-wall  flow,  especially  in  rough-wall  flow 
and  that  the  stress  due  to  small-scale  vortices  is  strongly  enhanced  over  the 
rough-wall.  An  alternative  way  that  we  attempted  to  extract  the  govern- 
ing event  associated  with  the  generation  of  turbulent  shear  stress  is  based 
on  LSE  (Linear  Stochastic  Estimation,  Adrian  et  al. , 1988).  The  extracted 
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(a)  Smooth  wall  at  ReT  = 395  in  (b)  Rough  wall  in  y+  = 10  ~ 50 
y+  = 10  ~ 50 

Figure  9.  Contour  lines  of  statistical  property  of  vorticity  on  a — /?  plane.  Solid 
lines  : positive,  dotted  lines  : negative. 


structures  in  positioning  the  detector  in  the  plane  parallel  to  the  wall  at 
y+  = 25  are  shown  in  Fig.8  which  demonstrate  that  streamwise  elongated 
vortices  are  the  key  structures  and  that  the  key  structure  does  not  depend 
on  the  position  of  the  detector,  as  in  (b)  which  is  for  the  detector’s  position 
unconditioned.  It  is  also  found  that  the  extracted  vortices  are  more  tilted 
from  the  wall  in  the  rough-wall  case  than  in  smooth  wall  case,  though  the 
pictures  are  not  shown  here.  It  is  concluded  that  such  a modulation  of  at- 
titude angle  of  vortices  is  substantial  in  enhancing  the  turbulence  mixing 
of  near-wall  layer,  in  rough-wall  flow.  In  order  to  know  the  tilting  property 
of  the  vortices,  weighted  average  of  a quantity  <j>  on  the  attitude  angle  of 
vortex  line  a,  f3  is  defined  as, 


P+ A/3  ra+ Aa 


f rp+Ap  r 

<t>  = / 

JVJ0-AP  Ja-Aa 


0(a,  (3)dadf3dV 


where  a = tan  1 i+>y / V +w|),  (3  — tan  1(ujx/u>z),  Aa  = 2.5°  ! $A (3  = 
5° ! % 

Figures  9 which  shows  the  average  for  cf>  = Q gives  the  peak  of  Q in 
a- (3  map  at  around  f 3 ss  90°,  which  corresponds  to  streamwise  small-scale 
vortices.  In  rough  wall  flow,  the  peak  shifts  to  larger  a than  in  smooth-wall 
flow,  indicating  that  the  vortex  structures  inclines  more  from  the  wall  in 
average. 


3.4.  FORMATION  OF  SMALL-SCALE  COHERENT  VORTICES  OVER  THE 
ROUGHENED  WALL 

Though  small-scale  coherent  vortices  are  found  to  be  key  structures  of  tur- 
bulence mixing  over  roughened  wall,  the  well-known  coherent  vortices  found 
in  smooth  wall  flow  are  different  however,  from  those  of  rough-wall  flow. 
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Figure  10.  Evolution  of  a hairpin  vortex  initially  introduced  in  a laminar  flow 
having  velocity  distribution  of  a turbulent  flow,  T+  =0,  120,  240.  Light  color  : 
Q+  > 0.0001,  Dark  color  : u'  < - 2 

Visualization  of  the  small-scale  coherent  vortices  demonstrates  diversity  of 
vortices.  In  the  region  over  the  top  of  roughness  element,  it  seems  that 
vortex  generation  is  strongly  associated  with  the  flow  separation  at  the 
upstream  corners  of  roughness  elements,  via  a shedding  of  strong  vortex 
sheets.  However,  complexity  of  flow  structure  does  not  allow  us  to  trace 
the  evolution  of  vortices.  Therefore,  to  investigate  the  fundamental  process 
of  vortex  formation,  an  evolution  of  a single  hairpin  vortex  is  traced  by 
conducting  numerical  simulation  in  which  a laminar  flow  having  mean  flow 
distribution  of  turbulent  flow  is  assumed  and  small  eddy  is  introduced  in  it 
at  an  initial  time.  The  introducing  baby  eddy  can  be  extracted  by  LSE  as 
proposed  by  Zhou  et  a/.  (1996).  Figures  10  show  the  evolution  for  the  non- 
dimensional  time  T+(=  Tu\/v 2)  = 0,120,240  and  depict  iso-surfaces  of 
Q+  and  low-speed  streaks.  The  vortex  introduced  initially  is  a pair  of  small 
quasi-streamwise  vortices.  At  the  beginning,  the  quasi-streamwise  vortices 
grows  quickly  to  a hairpin  vortex  and  extends  to  further  away  from  the  wall. 
At  T+  = 120,  the  secondary  hairpin  vortex  is  born  newly  at  the  upstream 
of  primary  hairpin  vortex,  i.e.,  regeneration  is  confirmed.  At  T+  = 240, 
the  group  of  vortices  becomes  matured  and  the  size  of  the  originally  intro- 
duced vortex  becomes  larger  than  half  width  of  the  channel.  Qualitatively, 
the  evolution  of  vortices  in  the  upper  layer  from  the  wall  is  similar  to  that 
of  smooth  wall  case  simulated  in  the  previous  work  (Zhou  et  al.  1996, 1999 
; Miyake  et  al.,  2000).  However,  in  the  layer  closest  to  the  wall,  a distinct 
pattern  of  vortex  generation  appears.  In  the  smooth  wall  case,  though  not 
shown  here,  a lot  of  quasi-streamwise  vortex  are  generated  over  the  flat 
wall,  most  probably  forming  a vortex  group.  In  rough  wall  case,  instead  of 
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such  a group  of  quasi-streamwise  vortices  around  the  legs  of  primaly  hair- 
pin vortex,  vortex  shedding  at  the  upstream  corner  of  roughness  elements 
occurs  induced  by  the  vortical  flow  from  upper  layer.  Although  the  ordi- 
nary regeneration  process  of  smooth  wall  flow  are  destroyed  in  rough  wall 
flow,  strong  vorticity  from  not  random  but  fixed  position  associated  with 
roughness  supplements  the  loss  and  maintain  turbulence.  This  character- 
izes the  rough-wall  turbulence  as  an  alternative  way  of  generating  vortices, 
suggesting  that  regeneration  of  quasi-streamwise  vortices  by  themselves  is 
not  always  unique  process. 

4.  Conclusions 

1.  The  turbulence  property  normally  found  in  the  layer  away  from  the  wall 
prevails  close  to  the  wall  and  large-scale  streaks  comes  closer  to  the  wall. 
It  is  suggested  that  a large-scale  mixing  takes  place  over  the  roughness. 

2.  The  major  effect  of  roughness  element  is  to  modify  the  mechanism  of 
turbulent  mixing.  Modulation  of  attitude  angle  of  vortex  which  is  tilted 
upward  from  the  wall  is  the  key  factor  in  enhancing  mixing  in  near-wall 
layer  over  the  roughness. 

3.  On  the  observation  of  simplified  simulation  and  the  visualized  instanta- 
neous flow  field  of  full  simulation,  further  away  from  the  wall,  the  mech- 
anism of  generation  of  coherent  vortices  is  similar  to  that  of  smooth  wall 
case.  Closer  to  the  top  of  roughness,  however,  quasi-streamwise  vortices  are 
not  observed.  The  vortex  sheet  shedding  at  the  upstream  corner  of  rough- 
ness causes  the  vortex  generation.  Thus  in  rough-wall  flow,  regeneration 
process  of  the  quasi-streamwise  vortices  by  themselves  is  suppressed. 
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